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<—| ' Abstract 

^ r~| Generalization of the Kac integral and Kac method for paths measure based on the 

-i— > . 

Levy distribution has been used to derive fractional diffusion equation. Application to 
nonlinear fractional Ginzburg-Landau equation is discussed. 

> 

1 Introduction 

■^j- ■ Kac integral [H El Ej appears as a path-wise presentation of Brownian motion and shortly 

^ ' becomes, with Feynman approach [4], a powerful tool to study different processes described by 

O 

the wave-type or diffusion-type equations. In the basic papers [HH], the paths distribution was 
based on averaging over the Wiener measure. It is worthwhile to mention the Kac comment 

5— i ' 

that the Wiener measure can be replaced by the Levy distribution that has infinite second 
and higher moments. There exists a fairly rich literature related to functional integrals with 
generalization of the Wiener measure (see for example (5j Ej). Recently the Levy measure 
was applied to derive a fractional generalization of the Schrodinger equation [7J [8] using the 
Feynman-type approach and expressing the Levy measure through the Fox function [9] 

In this paper, we derive the fractional generalization of the diffusion equation (FDE) from 
the path integral over the Levy measure using the integral equation approach of Kac. 
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2 Levy distribution 

Let us consider the transition probability P(x, t\x', t') that describes the evolution of the prob- 
ability density p(x, t) by the equation 

r+oa 

p(x,t)= / dx' P(x,t\x',t')p(x',t'), (1) 



where 

/+oo 
dxp{x,t) = l. (2) 
-oo 

The function P(x,t\x' ,t') can be considered as conditional distribution function. Then the 
normalization condition 

/+oo 
dxP(x,t\x',t') = 1 (3) 
-oo 

holds. Assume that P(x,t\x' ,t') satisfies the Markovian (semigroup) condition 

r+oo 

P{x, t\xQ, t ) = / dx'P(x,t\x',t')P(x',t'\xo,t Q ) (4) 



known also as the Chapman- Kolmogorov equation. 

In physical theories, the stability of a family of probability distributions is an important 
property which basically states that if one has a number of random variables that belong 
to some family, any linear combination of these variables will also be in this family. The 
importance of a stable family of probability distributions is that they serve as "attractors" for 
linear combinations of non-stable random variables. The most noted examples are the normal 
Gaussian distributions, which form one family of stable distributions. By the classical central 
limit theorem the linear sum of a set of random variables, each with a finite variance, tends to 
the normal distribution as the number of variables increases. All continuous stable distributions 
can be specified by the proper choice of parameters in the Levy skew alpha-stable distribution 
[TO] that is defined by 

r+oo 

L(x,y,a,/3,c) = — dpe~ ipx U(p,y,a, 0,c), (5) 
2vr y_ 00 

where 

U(p,y,a,(3,c) = exp(iyp- \cp\ a [l - i(3sign(p)$(a,p)}) , (6) 



and 

{tan(7ra/2), < a < 2, a ^ 1; 
(7) 
-(2/7r)logM, a = l. 

Here y is a shift parameter, /3 is a measure of asymmetry, with [3 = yielding a distribution 
symmetric about y. In Eq. ([6]), parameter c is a scale factor, which is a measure of the width 
of the distribution and a is the exponent or index of the distribution. 

Consider P(y,t'\x,t) as a symmetric homogeneous Levy alpha-stable distribution 

P(y,t'\x,t) = K(y-x,t'-t) = ^-J dpexp(ip(y -x) - (f -t)C a \p\°) , (0 < a < 2). (8) 

For a = 2, Eq. (jHJ) gives the Gauss distribution 

P(y,t / |x,t)= . 1 = expf- — — - t -(y-x) 2 ). (9) 

V ' ^ v/47rC 2 (t' - t) V 4C 2 (t'-t) Vy V V ; 

Eq. (jSJ) gives the function 

K( x ,t) = —[ dpexp(ipx-tC a \p\ a ) (10) 
2vr V / 

that can be presented clS cl Fourier transform 

Jjf( a: ,f)=^ 1 ( c ^«W a ) J (11) 



where 



For a = 2, Eq. (TTTjl gives 



^(/(P)) = TT" / dpe^f(p). (12) 



2tt 



if(a?,t) = -=i==expf--^-). (13) 



In the general case, the function K(x,t), given by Eq. (TTIi) . can be expressed in terms of the 
Fox ^-function El El HH G21 UHl 031 (see Appendix). 
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3 Fractional Kac path integral 



Let us denote by C[t a , tb] the set of trajectories starting at the point x a = x(t a ) at the time t a 
and having the endpoint Xb = x{tb) at the time tb- 
The Kac functional integral [21 [3l [15] is 

ftb 

V w x(t) exp( — / 

*C{t a ,t b ] 



[ V w x(t) expf- dTV(x(r))), (14) 



W(x b ,t b \x a ,t a ) -- 
where V(x) is some function, and 

V w x = lim \\ K{Ax kl At k )dx k . (15) 

For ffl3l) . expression fTl5|) gives 



n— >oo ■ 

fc=l 



n 



2V* = hmTT - dXk expf-j^L), (16) 

which is the Wiener measure of functional integration [15]. The integral ffT4l) is also called the 
Feynman-Kac integral. Using ffTUl) for a = 2, the path integral (|14l) can be written as 

If n 
W(x b ,tb\x a ,t a ) = lim / dxidp 1 ...dx n dp n exp\^(ip k Ax k - At k [C 2 pl + V(x k 

n-oo (27T) n y M 2n f-'V 



(17) 



where the time interval [t a , tb] is partitioned as 



t k — t a + k— -, to — ta, t n — t b , (18) 

n 

and 

Ax k = x k+1 -x k , At k = t k+1 -t k , x k = x(t k ), p k =p(t k ). (19) 
The functional integral (I17p can be rewritten as 

W(x b , t b \x a , t a ) = ! VxVp exp( j dt ipx - C a p 2 - V(x) ) . (20) 



where 

n U d 

Vx = lim | | dx k , Up = lim I I . (21) 

n. — +nn n — >r>r, 



k=l k=l 



The Kac functional integral in the form ( 1201) is a classical analog of the Feynman phase-space 
path integral, which is also called the path integral in Hamiltonian form. 



For the fractional generalization of Wiener measure (JT5l) and Kac integral (1T4"1) . we consider 
K(x,t) given by §10f). Substitution of §10§ into 

„ n 

W(x h t b \x a ,t a ) = lim / ]ldx k K(Ax k ,At k )exp(-At k V(x k )) (22) 
n ^°° Jul* ^ v J 



with 



K(Ax k , At k ) = ^- J dp k exp(ip k Ax k - At k C a \p k \ a ^J, (0 < a < 2), 



gives 



W(x b ,t b \x a ,t a 



lim 

n— >oo 



n d d n 

II — \ — - exp^2(ip k Ax k - At k [C a \p k \ a + V(x k )] 



k=l k=0 

Similarly to (|20|) . (I2T1) this expression can be written as 



W(xb, h\x a , t a ) = / VxVp exp^/ <it ipx — C Q |p| Q — V(o;) J. 



This expression is a fractional generalization of (1201) . 
If we introduce formally imaginary time such that 

dx dx 



ix = i- 



dt ds 

then (1251) transforms into the Feynman path integral with a generalized action [7J Ej 



dt 



px-C a \p\ a -V{x) 



as an action. Hamiltonian-type formal equations of motion are 

— = N a \p\ a - 1 dp - 

ds " ' ds dx ' 



(23) 



(24) 



(25) 



(26) 



where N a = aC a sign{p). 
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4 Fractional diffusion equations 

It is known that the Kac integral (1141) can be considered as a solution of the diffusion equation 
[21 IS] . Let us derive the corresponding diffusion equation for the fractional generalization of 
the Kac integral (|25|) . 

In (j25p the integration is performed over a set C[t a ,tb] of trajectories that start at point 
x a = x{t a ) at time t a and end at point Xb = x(tb) at time U- For simplification, t a = 0, x a = 0, 
and tb = t, Xb = x are used. In particular, we can consider two following cases of C[t a ,tb\. 

(1) The set C/[0,t] consists of paths for which both the initial and final points are fixed. 
The integration over this set obviously gives the transition probability 



The conditional fractional Wiener measure corresponds to the integration over the set C/[0,£] 
of paths with fixed endpoints: x a = 0, 

(2) If we consider a set C a [0,t] of trajectories with arbitrary endpoint Xb = x, the measure 
is called the unconditional fractional Wiener measure. This measure satisfies the normalization 
condition 




V w x = K(x b - x a , t b - t a ) = P(x b , t b \x a , t a ), 



or 





(27) 



since it is a probability that the system ends up anywhere. 
For simplification, we introduce the notation 




(28) 



and define the field 



u(x,t) = W(x,t\0,0). 



(29) 



For the fractional Kac functional integral, we have with respect to (|27|) . 




(30) 



6 



Using notations (I28j) . (12911 . expression ( 1251) for £ a = 0, x a = 0, and U = t, Xb = x can be 
presented as 



t) 



/ XV^expf — / drV(x(r))J = / T> w xZ[x,t]. 
Jc f [o,t] ^ Jo ' Jc[o,t] 



(31) 



To derive a fractional diffusion equation, we use the identity [T5 



cxp 



dr^(x(r))) = 1 - / dr U(x(r))exp 



dsV(x(s] 



(32) 



Equation ( |32jl can be proved by using differentiation by t, and the value of the constant is found 
from the condition of coincidence of both sides for t — 0. For the notation (I29p . identity 
has the form 



Z[x,t] = 1 - / dr V(x(t))Z[x,t] 



(33) 



Equation (|33|) can be integrated with respect to the conditional fractional Wiener measure: 



V w x Z[x,t] 



c f [o,t] 



V w x 1 



C f [0,t] 



T>wx I dr 



c f [o,t] 



V(x(t))Z[x,t] 



(34) 



Changing the order of the integration in the second term in the right hand-side of (|34|) . we get 



V w x / dr 



C f [0,t] 



V(x(t))Z[x,t] 



dr 



V w x 



dr 



dx T 



V w x 



dr 



JCf 



>Cf[0,T] 

dx T V(x(t)) 



c f [o,t] 
V w x 



V(x(t))Z[x,t] 



C f [r,t] 



V{x{t))Z[x,t 
/ V w xZ[x, r] / V w x. 

JCf[0,T] JC f [r,t] 

The first term in the right hand-side of ( 1341) gives 

V w xl= lim / TT dxk K(Axk, Atk) — 

„ n 

= lim / TT dajjfc (Aa; fe) At fe ) = i^(x 6 -x a ,t b — t a ). 
Using (1291 . (|36|) . and (1351) . Eq. (|34|) gives the integral equation 



(35) 



+oo 



u(x,t) = K(x,t) — dr dx T V(x T ) u(x T , r)K(x — x T , t — r). 



(36) 



(37) 



For this equation there exists the infinitesimal operator C a (generator) of time shift such that 

du(x, t) 



dt 



C a u(x,t). (38) 



Using (EH) and (ITTj) . we obtain 



d a 1 /"* 

C a u(x, t) = C a — — j— u(x, t) — lim - dr I dyK(x — y,t — T)V(y)u(y,r), (39) 



d\x\ a v ' ' t-at 

where d a /d\x\ a is a fractional Riesz derivative [T6 l fT7 l fT8 l fT9] of order < a < 2 that is defined 
by its Fourier transform 



d° 



d\x\ a 

where 



u(x, t) = T- 1 (\p\ a u(P, t ))=^J d P \P\^ *) e ~ iPX ' ( 4 °) 



-+oo 

u(p,t) = I dxu(x,t)e ipx . (41) 



The initial condition K(x, 0) = 5(x) gives [15 



1 



lim- / dr / dyK(x-y,t-T)V(y)u(y,r) = V(x)u(x,t). (42) 
Then ( 1391) gives 

^ = ^^"^)- (43) 
This generator is an operator of fractional differentiation of order a. 
As a result, we obtain 

du(x,t) d a u(x,t) Ir , , , , .... 

-k^ =c °^^- niMa; - t) ' (44) 

which is a diffusion equation with fractional coordinate derivatives. For a = 2, Eq. ( 1441) is the 
usual diffusion equation. 

It is worthwhile to mention that the way of obtaining fractional equation (|44|) is based on 
the exploiting the properties of integral equation ( 1371) . while the expansion of exponents in (j21 
over small Atk has been used in [8] for Feynman path integral. 
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5 Fractional diffusion equations by Kac approach 

It is useful also to derive the fractional diffusion equation from (TJ3J) using Kac approach de- 
scribed in Sec. 4. of [2]. 

The mathematical expectation value of Z[x,t] is defined as 



£^exp(- jf drV(x(T)))^ = ^ £> w x ex p(^ drF(ar(r)). 



Using the expansion 



we get 



QOM) = V / V w x expf / rfrU( 



x(t)) 



(45) 



expf- fdrV{x{r))) = f) /' drV(x(r))) * (46) 

Wo / ^ m! vy y 



£?/exp(- f drV{x{r)))\ = ^ LiT / V w x [f drV{x{r)) m . (47) 

The expression (H7I) can be presented as 

^/expf- f drV(x(r)))\ = V(-l) m dxQ m (x,t), (48) 
\ Wo y / m =o •'-«» 

where 

Q m (x,t) = — f V w x( [ drV(x(r))Y. (49) 
™ yo/to,*] Vo ' 

These functions (1491) satisfy the recurrence equations [2] 

/•t r+oo 

Q m +i(x,t) = dr dyK(x -y,t-r)V(y)Q m (y,T), (50) 

JO J-oo 

and 

QoOM) =A"(x,t). (51) 

Let us introduce 

oo 

Q{x,t) = ^2{-l) m Q m (x,t). (52) 

m=0 

Then 



T>wx expl / drV(x(r) 
>[o,t] vo 



(53) 



and 



s/exp(-^ dry(x(r)))\ = y dxQ(x,t). (54) 
It follows from (1501) and (i5"T|) that the field Q(x,t) satisfies the integral equation 

Q(x,t) = Q (x,t) - dr dyK(x-y,t-T)V(y)Q(y,r). (55) 

</0 </-oo 

There exists an infinitesimal operator C a of time shift such that 

^^- = £. a Q(x,t). (56) 

Using (1501) dlSD, and (TTTi) . this generator can be expressed through a fractional differential 
operator 

£ a Q(x, t) = C a — — Q(x, *) - lim - / dr / dyAT(x - y, t - r)V(y)Q(y, r). (57) 
The initial condition K(x, 0) = 5(x) gives similar to ( |42l) 

lim - / dr / eft/i^z - y, t - r)V(y)Q(y, r) = V(x)Q(x, t). (58) 



As a result, we obtain 



dQ{x,t) _ c PQ(x,t) 



dt d\x 

which is fractional diffusion equation that coincides with (jHJ). Then 



Q(x,t) =W(x,t\0,0) = I V w xexp(-[ drV(x(r))). (60) 

JctlO.t] ^ Jo ' 



'Cf[0,t] v Jo 

Using fl52j), the approximate solution of (JHJ) can be presented as 

u(x, t) ps Qo(^7 - Qi(x, t) + Q 2 (x, t) 

-t /"+oo 



ft r+oo 

= K(x,t)- dr dyK{x-y,t-r)V{y)K{y,r)+ 

JO J-oo 

+ dr dt' dy dy'K(x-y,t-T)V(y)K(y-y',T-t')V(y')K(y',t>). (61) 

JO JO J-oo J-oo 



for small enough V(x). 
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6 Nonlinear fractional equations 



Equations (jSJ) and (ISTJj) are linear equations with respect to the fields u(x, t) and Q(x,t). 
In general, nonlinear equations can be derived from the functional integral over the space of 
branching paths (see [21] and Sec. VI. 4. of |20j). Note that Feynman path integral over the 
branching paths has been suggested in [22] (see also [2"3"] |2"4] ) . The multiplicative representations 
of nonlinear diffusion equations are also considered in [25] [26], [27j . As an example of nonlinear 
diffusion equation, which can be derived from integrals over the branching paths, is an equation 
with the polynomial nonlinearity [20, 21 J : 

m 

U(u) = J2^kHx,t)] k . (62) 

k=2 

Using fractional Kac integral over the branching Levy paths [28, 29J, a nonlinear generalization 
of fractional equation (jUJ) can be derived in the form 

= C -^W^ ~ V V)<*> t) + E *)]*■ ( 63 ) 

I 3 '' k=2 

For example, fractional equations with cubical nonlinearity can be obtained 

^r 1 = ^^F 1 ~ v{x)u{x ' t] + asMx ' t)? - (64) 

Equation ( 1641) is the fractional generalization of the Gross-Pitaevskii equation [30] [31]. For 
V(x) = const, Eq. (IM1) is fractional Ginzburg-Landau equation that is suggested in [32] (see 
also [33] [31]) to describe complex media with fractional dispersion law. 
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Appendix: Fox function representation for K(x, t) 

In this section, we use the results of the paper [7j (see also [8]) to demonstrate how the function 
K(x, t) defined by Eq. fflQ]) can be expressed in the terms of the Fox if-function [9] [TTJ [121 
[13j [T3]. The Fox function representation of K(x,t) can be considered as a fractional analog 
of expression ([TBI . To present K(x,t) in terms of the Fox if-function, we consider the Mellin 
transform of fflQ]) . Comparing of the inverse Mellin transform with the definition of the Fox 
function , we obtain an expression in terms of Fox //-function. 

Using the relation K(x,t) = K(—x,t), it is sufficient to consider K(x,t) for x > only. 
The Mellin transformation of fflOl) is 

A f°° 1 f°° f + °° / \ 

K (s, t) = / dxx s K(x,t) = — / dxx dp exp[ipx - C a \p\ a t) . (65) 

Jo 2tt j y_ oc V / 

Changing the variables 

p ^ (Cat)- 1 '" Tl, X^{C a tf' a t 

A 

we present K (s,t) as 

K (s,t) = — {(C a t)^) s - 1 / d^*- 1 / dr)^-H\ ( 66 ) 

The integrals over g?£ and c??7 can be evaluated by using the equation [13]: 

fVr- 1 r^e^ = ^-sin^f^r( S )r(l-^), (67) 
7o Jo 8-1 2 \ a / 

where s — 1 < a < 2 and T(s) is the Gamma function. 

Inserting of flBTj) into f[66l) and using the relations 

T(l - ^) = -zr(-z), r(*)r(l -z)=n/ sin ttz, (68) 

we find 

^ (M )4( (0 „ t) v. r iWrg.. (69) 

Then the inverse Mellin transform of ([69]) is 



c+ioo c+ioo _^ 
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where the integration contour is the straight line from c — too to c + too with < c < 1. 
Replacing s by — s, we get 

-c+ioo 

(71) 



K(x,t) = -(C a t)- 1 / a ^- 
a lixi 



r(^)r(if)" 



The integration contour may be deformed into one running clockwise around [— c, oo). Com- 
parison with the definition of the Fox ff-function [91 [11] [12] gives 



1 

a 



K(x,t) = -(C a ty 1/a H 1 / 2 



{C a t)- 1/a x 



;i-l/a,l/a),(l/2,l/2) 



(0,1), (1/2, 1/2) 

Using the properties of the Fox H- function [91 [Til EE] , we obtain 



(72) 



lf(a:,*) 



ax 



-if. 



i.i 

2.2 



(C Q t)- 1/a |a;| 



;i,l/a), (1,1/2) 
(1,1), (1,1/2) J 



(73) 



Let us show by analogy with [7] (see also [8]) that Eq. (1731 includes as a particular case at 
a = 2 the well known Gauss distribution (jTBl . Assuming a = 2 in Eq. (IT3"]) . 

;i, 1/2), (1,1/2) - 



Ji-(x,t)| a=2 = /f 2 1 ; 2 1 (C^)- 1 / 2 ^! 

The series expansion of the function f[T4"|) gives 

^,t)U= 2 = ^(c 2 t)- 1/2 ^(-(c 2 t)- 1 / 2 ) 



(1,1), (1,1/2) 



(74) 



k \X\ 



k=0 



k\ r(±^) 



(75) 



Substituting of k — > 2/ into (!75|) . and using 



rl .2"' 



the function K(x,t) can be rewritten as 



(2) 2/ Z! 



(-i)W 



(76) 



K(x,t)| 



o=2 



x 



ACot 



(77) 



Thus, it is shown that (|T3|) can be derived from equation (1731 with a = 2. 
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